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TURBULENT  HEAT  TRANSFER  IN 
LARGE  ASPECT  CHANNELS 


F.D.  Haynes  and  G.D.  Ashton 


INTRODUCTION 

Turbulent  heat  transfer  in  round  tubes  has  been 
studied  extensively  and  correlations  between  heat 
transfer  and  flow  rate  are  well  established. 
Turbulent  heat  transfer  in  rectangular  channel 
geometries  has  been  studied  comparatively  less. 

This  study  was  motivated  by  the  authors’  interest 
in  heat  transfer  from  a river  to  its  ice  cover;  the 
cross-section  geometry  of  the  river  is  better 
specified  by  an  aspect  ratio  B/D  where  B is  the 
channel  width  and  D is  the  flow  depth.  This 
report  presents  new  data  for  a closed  channel  flow 
with  BID  = 10  and  compares  these  data  with  the 
earlier  data  of  Ashton  (1971)  and  Hsu  (1973).  The 
data  from  the  present  study  agree  reasonably  well 
with  the  predicted  data  of  Petukhov  and  Popov 
(1963). 


EXPERIMENTS 

The  experiments  were  conducted  in  a rectangular 
channel  with  width  B = 0.254  m and  flow  depth  D 
= 0.0254  m.  A schematic  of  the  apparatus  is  shown 
in  Figure  1 and  a cross-section  schematic  in  Figure 
2.  The  lower  boundary  of  the  apparatus  was  steel 
and  enclosed  a flow  of  trichloroethylene  (TCE) 
circulated  below  the  water  channel  at  an  inlet 
temperature  of  about  -28  °C.  The  flow  depth  for 
the  refrigerant  was  0.0254  m.  The  top  and  sides 
of  the  water  channel  were  constructed  of  wood  and 
the  entire  channel  section  was  insulated  with  ure- 
thane and  fiberglass.  The  test  section  was  3.66  m 
long.  At  the  inlet,  a 4000-W  immersion  heater  was 
controlled  by  a thermostat  to  allow  control  of  the 
water  temperature.  The  bulk  temperature  of  the 
water  flow  was  from  about  3°C  to  7°C.  Water  flow 
rates  were  from  0.64  x 10'3  to  4.5  x 10'3  m3  s’1 
with  corresponding  mean  flow  velocities  from  0.1 
to  0.7  m s '.  Temperatures  of  the  water,  wall,  and 
refrigerant  were  measured  at  three  stations.  The 


distance  from  entry  to  the  first  station  was  1 .676  m. 
Based  on  the  results  of  Hatton  and  Quarmby  (1963), 
the  First  station  was  well  beyond  the  distance  at 
which  entry  effects  are  significant.  The  distance 
between  stations  was  0.9144  m.  After  flow  ad- 
justments were  made,  the  system  was  allowed  to 
stabilize  for  about  30  minutes  before  temperature 
data  were  collected. 

Four  thermistors,  each  calibrated  at  the  water- 
ice  triple  point,  were  positioned  (see  Fig.  2)  at  each 
of  the  three  stations  and  the  thermistor  signals  were 
read  with  a Keithley  Model  172  Digital  Multimeter 
in  a forward-reverse-forward  sequence.  Based  on 
signal  ranges  used  and  manufacturer’s  specifications, 
absolute  temperatures  were  determined  to  about 
0.01°C  and  differential  resolutions  to  about 
0.005°C.  The  heat  transfer  rate  was  determined 
from  measurements  of  the  decrease  in  temperature 
of  the  flow  between  measuring  stations  with  the 
heat  flux  assumed  to  be  only  from  the  flow  to 
the  bottom  steel  plate.  This  last  assumption  was 
verified  by  order-of-magnitude  calculations.  The 
experiments  were  also  conducted  at  wall  tempera- 
tures above  0°C  so  that  no  ice  formed  on  the  wall. 


RESULTS 


A total  of  27  tests  were  conducted.  The  results 
are  summarized  in  Table  1 in  terms  of  the  Prandtl 
number,  Pr,  Reynolds  number,  Re,  and  Nusselt 
number,  Nu,  defined,  respectively,  by 


Ul\p 


Nu  = 

*(roc-r0) 


where  U is  the  mean  flow  velocity,  q is  the  heat  flux 
from  the  flow  to  the  wall,  and  To-T0  is  the  differ- 
ence between  the  bulk  water  temperature  To  and 
the  wall  temperature  T0.  p,  p,  k,  and  Cp  are  the 
dynamic  viscosity,  density,  therm*!  conductivity,  and 
specific  heat  capacity  of  the  water  and  were 
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Table  2.  Summary  of  experimental  results  of 
Ashton  (1971)  and  Hsu  (1973) 


Ashton  (1971) 


evaluated  at  the  bulk,  or  mixed  flow  temperature. 

The  hydraulic  diameter  £>h  is  defined  in  the  conven- 
tional way  by  Dh  = 4 A IP  where  A is  the  cross-section 
area  of  the  flow  and  P is  the  wetted  perimeter.  For 
the  present  experiments,  Dh  = 1.818  D.  All  of  the 
present  experiments  were  in  a narrow  range  of  Prandtl 
number  9.90  < Pr  < 12.28  with  a representative  Pr 
for  most  of  the  data  of  Pr  =“  1 1 .8 

Ashton  (1971)  and  later  Hsu  (1973)  investigated 
turbulent  heat  transfer  at  the  bottom  of  an  open 
channel  flow  in  a flume  0.61  m wide  and  12.2  m 
long;  most  experiments  were  conducted  with  a 0.1  S-m 
flow  depth.  They  froze  an  ice  slab  onto  the  bottom 
of  the  flume  and  measured  the  heat  flux  by  cal- 
culating the  change  in  thickness  of  the  ice  at  the 
center  of  the  channel  and  well  downstream  from  the 
inlet.  Their  test  results  are  given  in  Table  2.  In  their 
experiments,  Dh  = 2.667  D0  where  D0  is  the  flow 
depth.  In  their  experiments,  the  range  of  Prandtl 


number  was  quite  small  (1 2.34  < Pr  <1 1 3.62)  with 
a representative  value  of  Pr  =“  1 3.3. 

DATA  INTERPRETATION 

All  the  data  are  plotted  in  Figure  3 as  a function  of 
Re  and  Nu.  Three  of  Hsu’s  data  points  overlap  the 
data  of  the  present  study.  For  comparison,  two 
analytical  predictions  are  also  plotted  in  Figure  3. 

The  formula  of  Petukhov  and  Popov  (1963)  corre- 
lates turbulent  heat  transfer  in  round  tubes  quite 
well  (t6%)  for  a wide  range  of  Re  and  Pr  [see,  e.g., 
Karlekar  and  Desmond  (1977,  p.  351  )1  ; it  underpre- 
dicts the  Hsu  (1973)  and  Ashton  (1971)  data  by  about 
1 5%,  and  the  data  of  the  present  study  by  about  35%. 
The  formula  of  Petukhov  and  Popov  (1963)  is 
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Nu  = . 


(// 8)  Re  Pr 


Figure  3.  Heat  transfer  as  a function  of  Reynolds  number. 

-0.015  for  the  Prandtl  numbers  of  the  present  data 
sets,  so  that  the  Shibani-Ozisik  formula  becomes 


1.07+12.7  (//8),/2  (Pr2/3-1) 


(D 


where  / has  been  calculated  using  the  Filonenko 
equation  [see  Karlekar  and  Desmond  (1977,  p.  352)] 


1.82  log]0  Re-1.64]- 


(2) 


We  also  note  that  the  PetuknoWteaQv  formula  yields 
Nusselt  numbers  somewhat  greater  (~2ft%)»thaathe 
commonly  used  Dittus-Boelter  equation  or  the  CoP* 
burn  equation  in  the  range  of  Reynolds  numbers  of 
the  Hsu  and  Ashton  data  sets. 

Shibani  and  Ozisik  (1977)  presented  an  analytical 
solution  for  turbulent  heat  transfer  between  parallel 
plates  with  top  and  bottom  cooling  and  correlated 
their  analytical  predictions  in  the  range  0.1  < Pr  < 
104  and  1 04  < Re  < 1 06  by 


Nu  = 12+0.03  Re*1  Pr*2 


(3) 


Nu  = 12+0.03  Re0-865  Pr0-33 


(6) 


Equation  6 is  also  presented  in  Figure  3 and  is  seen 
to  overpredict  the  Nusselt  number  by  about  50%. 

Several  reasons  are  offered  for  the  differences 
between  the  data  and  the  predictions.  First,  it  is 
difficult  to  achieve  the  idealized  flow  conditions 
implicit  in  the  analyses.  In  particular,  we  suspect 
the  steel  plate  had  some  roughness  in  the  present 
genments  and  in  the  Reynolds  range  of  these 
expenriTenlfcJtmcould  well  have  contributed  to  a 
higher  Nusselt  nufiVt>w*^$econd,  the  aspect  ratio  fi/Oh 
varied  from  1.5  for  the  AshtoiWsu.4gta  sets  to  5.4 
for  the  present  study,  while  the  Petukhov^Pbp#»^for- 
mula  is  applicable  to  a round  cross  section  roughly 
equivalent  to  a B/Dh  = 1 and  the  Shibani-Ozisik 
formula  is  applicable  to  a B/Dh  = “>. 


where 
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a1  - 0.88  • 


0.24 

(3.6+Pr) 


(4) 


a2  = 0.33+0.5  e °-6Pr 


(5) 


For  the  Prandtl  number*  of  the  present  data  sets, 
the  second  term  in  eq  5 is  negligible  compared  with 
the  first.  In  eq  4,  the  second  term  is  approximately 
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